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Abstract 



We investigate rest-frame near-infrared (NIR) morphologies of a sample of 139 galaxies with M s > 
1 x 10 10 M at z=0.8-1.2 in the GOODS-North field using our deep NIR imaging data (MOIRCS Deep 
Survey, MODS). We focus on Luminous Infrared Galaxies (LIRGs), which dominate high star formation 
rate (SFR) density at z <~ 1, in the sample identified by cross-correlating with the Spitzer /MIPS 24/im 
source catalog. We perform two-dimensional light profile fitting of the z ~ 1 galaxies in the -fsT s -band (rest- 
frame J- band) with a single component Sersic model. We find that at z ~ 1, ^90% of LIRGs have low 
Sersic indices (n < 2.5, similar to disk-like galaxies) in the its-band, and those disk-like LIRGs consist of 
t— » . ^60% of the whole disk-like sample above M s > 3 x 10 10 Mq. The z ~ 1 disk-like LIRGs are comparable 

' or ~ 20% small at a maximum in size compared to local disk-like galaxies in the same stellar mass range. 

If we examine rest- frame UV-optical morphologies using the HST / ACS images, the rest-frame S-band 
sizes of the z ~ 1 disk-like galaxies are comparable to those of the local disk-like galaxies as reported by 
previous studies on size evolution of disk-like galaxies in the rest-frame optical band. Measuring color 
gradients (galaxy sizes as a function of wavelength) of the z ~ 1 and local disk-like galaxies, we find that 
the z ~ 1 disk-like galaxies have 3-5 times steeper color gradient than the local ones. Our results indicate 
(f) . that (i) more than a half of relatively massive disk-like galaxies at z ~ 1 are in violent star formation epochs 

observed as LIRGs, and also (ii) most of those LIRGs are constructing their fundamental disk structure 
, vigorously. The high SFR density in the universe at z ~ 1 may be dominated by such star formation in 

' disk region in massive galaxies. 
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1. Introduction Sanders & Mirabel 1996 for a review on LIRGs). Their 
, luminous IR emission is thermal dust emission in the mid- 
*_h ' It has been well known that the star formation rate IR (MIR) to far-IR (FIR) wavelength caused by absorbing 
(SFR) density of the universe (cosmic SFR density) in- UV photons from intensive star formation and/or lumi- 
creases by an order of magnitude from the present to z ~ 1 nous active galactic nuclei (AGNs). IR luminosity of 10 11 
(Lilly et al. 1995; Madau et al. 1996; Hogg et al. 1998; L Q corresponds to the SFR of ~ 15 M Q yr" 1 (Kennicutt 
Flores et al. 1999; Haarsma et al. 2000; Hopkins 2004; 1998), which is several to ten times larger than that of lo- 
Hopkins & Bcacom 2006). What kind of galaxies are con- cal normal galaxies (Brinchmann et al. 2004). Although 
tributing to such a large cosmic SFR density at z ~ 1 is in the local universe LIRGs are rare and have only ~ 5% 
still an open question. contribution to the IR luminosity density, the number den- 
In recent years with the Spitzer Space Telescope, many sity of LIRGs increases with redshift, and 70% of the IR 
studies focusing on understanding the properties of dis- luminosity density is in LIRGs at z ~ 0.7 (Le Floc'h et al. 
tant star- forming galaxies have been conducted. The cos- 2005). It is important to understand triggering processes 
mic SFR densities at high redshifts are dominated by in- of the LIRGs to unveil the physical cause of the evolution 
frarcd bright galaxies, especially Luminous Infrared (IR) of the cosmic SFR densities. 

Galaxies (LIRGs) which have IR (8-1000 fxm) luminosi- In the local universe, most of LIRGs show irregular mor- 

tics, Lir, of 10 11 -10 12 Lq (Flores et al. 1999; Takcuchi phologies, indicating galaxy-galaxy interaction (Sanders 

et al. 2005; Le Floc'h et al. 2005; Caputi et al. 2007; & Mirabel 1996; Sanders & Ishida 2004). Even among 

Perez- Gonzalez et al. 2008a; Magnelli et al. 2009; see the LIRGs which have spiral morphology apparently, a 
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significant fraction of them shows bar structure (Wang 
et al. 2006). Such morphological properties indicate that 
gas is pushed into the nuclear region and nuclear starburst 
is triggered. 

On the contrary, at high redshifts (0.7 < z < 1), the 
Hubble Space Telescope (HST) imagings show that most 
of LIRGs have spiral (or late-type) morphology without 
any clear sign of merging and/or interaction (Zheng et al. 
2004; Bell et al. 2005; Melbourne et al. 2005; Lotz et al. 
2008). Those galaxies also have no bar structure (Zheng 
et al. 2005). Moreover, in the last couple of years, some 
studies dedicated on MIR spectra of such distant LIRGs 
found their spectral shape indicates cool dust tempera- 
ture similar to local spiral galaxies with lower IR lumi- 
nosity (L m < 10 10 - 5 Lq) (Zheng et al. 2007; Symeonidis 
et al. 2009; Seymour et al. 2010). More recently, Takagi 
et al. (2010) reported that most of polycyclic aromatic 
hydrocarbon (PAH)-selected galaxies with Lir > 10 11 L Q 
at z ~ 1 show PAH-to-total IR luminosity ratio similar to 
that of less luminous starburst galaxies using an AKARI 
multi-wavelength MIR photometry. Those distant LIRGs 
show similar properties to local massive spiral galaxies ex- 
cept for their high IR luminosities. 

Previous morphological studies of LIRGs at high red- 
shifts were conducted in the rest-frame optical range. 
However, the high IR luminosity of LIRGs indicates the 
presence of a large amount of dust, and then morpho- 
logical analysis in Near-IR (NIR) range less affected by 
dust-extinction is necessary to examine existence of dust- 
obscured structures such as nuclear starburst seen in local 
LIRGs or bars invisible in optical. NIR wavelength is also 
suitable to reveal the distribution of old stellar popula- 
tions which dominate stellar mass of galaxies. As a pio- 
neering study in NIR, Melbourne et al. (2008) carried 
out targeted observations of 15 LIRGs at 0.4 < z < 1.2 
with Adaptive Optics (AO) in the if-band, and evaluated 
their morphologies visually. It was found that two thirds 
of the LIRGs look disky apparently while major mergers 
are ~ 30% at a maximum. They also found that only 
one LIRG has a clear bar structure. The distant LIRGs 
are different from local LIRGs in morphology, and star 
formation in the distant LIRGs seems to be smoothly dis- 
tributed in their disk. 

However, the morphological study by Melbourne et al. 
(2008) is based on the visual classification. Further NIR 
morphological studies with quantitative evaluation are 
necessary. In addition, although they observed normal 
galaxies (non-LIRGs) at the similar redshifts, it is still 
crucial to locate LIRGs in field galaxies in the same red- 
shift range. 

In this paper, we construct a stellar mass- and volume- 
limited, NIR morphological catalog of galaxies at z ~ 1 
independent of IR luminosity using deep iT s -band data, 
and locate LIRGs among z ~ 1 galaxies. We evaluate the 
galaxy morphology quantitatively using two-dimensional 
light profile in both K s - (rest-frame J-) and optical (rest- 
frame U- and V-) bands. Comparing the rest-frame NIR 
and UV-to-optical morphologies, we discuss distributions 
of old stellar population and currently star-forming re- 



gion in the high rcdshift LIRGs and non-LIRG galaxies. 
Throughout this paper, we use a A-CDM cosmology with 
Q m = 0.3, fl A = 0.7, and H = 70 km s" 1 Mpc" 1 . In 
this cosmology, 1 arcsec corresponds to 8.01 kpc at z = 1. 
We refer EST j Advanced Camera for Surveys (ACS) fil- 
ters F435W, F606W, F775W, and F850LP as B435, V 606 , 
1775, and zsso-bands, respectively. 

2. Sample of z ~ 1 Galaxies 

2.1. K s - selected stellar mass catalog from MODS Deep 
Data 

We construct a stellar mass-limited z ~ 1 galaxy sam- 
ple using the deep catalog of Kajisawa et al. (2009, here- 
after K09). The catalog contains 3203 if s (Vega) < 24 
mag galaxies in 28 arcmin 2 from MOIRCS (Multi-Object 
InfraRed Camera and Spectrograph, Suzuki et al. 2008) 
Deep Survey (MODS) ultra-deep image obtained with 
MOIRCS attached to the Subaru Telescope (Kajisawa 
et al. 2006). It covers a part of the Great Observatories 
Origins Deep Survey North (GOODS-N, Giavalisco et al. 
2004) field including Hubble Deep Field North region 
(HDF-N) . Utilizing publicly available multi-wavelength 
(U B i35 V 60 6 1775 z S50 J H K s , 3.6/xm, 4.5/xm, and 5.8/zm) 
data and spectroscopic rcdshift (z spec ) catalogs in the 
GOODS-N field, K09 constructed a catalog of rcdshift and 
stellar mass for the -Kg-selected sample. The redshifts and 
stellar masses were estimated by fitting model spectral 
energy distributions (SEDs) obtained from several popu- 
lation synthesis models to the SEDs of the sample. The 
Salpetcr Initial Mass Function (IMF) in stellar mass range 
of 0.1 - 100 M Q (Salpetcr 1955) is assumed in the SED 
fittings. In this paper, we adopt the photometric redshifts 
(zphot) and stellar masses (M s ) obtained with Galaxev 
population synthesis model (Bruzual & Chariot 2003). 
The comparison between the photometric and spectro- 
scopic redshifts of the galaxies which have the spectro- 
scopic data shows that the photometric redshift estimates 
have a good accuracy with median and standard deviation 
of 5z = (zphot - z sp cc)/ (1 + z spcc ) of -0.002 and 0.072 and 
the fraction of the catastrophic failure (8z > 0.5) of 4.2% 
for galaxies with z S p OC = 0.8-1.2. Note that the photomet- 
ric redshifts are used only for galaxies with no spectro- 
scopic identifications in order to maximize the reliability 
of the redshift catalog. The stellar masses were calcu- 
lated using the best- fit stellar mass-to-luminosity (M s /L) 
ratio in the observed if s -band. Typical uncertainties of 
the stellar masses estimated from the probability distri- 
butions in the SED fitting (Figure 2 in K09) show that 
the stellar mass errors at z < 1.5 are less than 0.1 and 
0.05 dex for stellar masses M s > 10 10 M Q and M s > 10 11 
Mq, respectively. Those errors include the uncertainty 
of the photometric redshift for galaxies with no spectro- 
scopic redshift. It should be mentioned that systematic 
and random uncertainties (typically 0.2 dex) associated 
with stellar population models and IMFs are not included 
in these values. 

Figure 1(a) shows the stellar masses of galaxies in the 
MODS deep region at z = 0.6-1.4. The MODS stellar 
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Fig. 1. Galaxies around z = 1 in the MODS deep region, (a) Stellar mass distribution of galaxies at z = 0.6-1.4. The dashed lines 
separate our z ~ 1 (z = 0.8-1.2) galaxies from the Ks-selected sample. The soKd line indicates the stellar mass limit corresponding 
to Ks(Vega) = 24 mag (Kajisawa et al. 2009, K09). The dotted line indicates a stellar mass limit for our morphological sample 
selection (M s = 1 X 10 1() Mq) defined in Section 4.1. Large filled circles (red) represent MIPS 24 ^im-detected (f24 > 80 i*Jy) 
sources. Chandra X-ray-detected sources are marked with large open circle, (b) Stellar masses and K B magnitudes of galaxies at 
z = 0.8-1.2. (c) Stellar masses and rest-frame U — V colors of galaxies at z = 0.8-1.2. In all panels, the uncertainties of the stellar 
mass estimates are less than 0.1 and 0.05 dex for galaxies with M B = 10 10 — 10 11 Mq and with M B > 10 11 Mq, respectively. 



mass catalog is complete down to at least M s ~ 10 9 Mq 
at these redshifts (shown with the solid line in the fig- 
ure). In this paper, we focus on galaxies at 0.8 <z < 1.2, 
which are plotted with black dots in the figure. For those 
z ~ 1 galaxies, we plot the distributions of the K s mag- 
nitudes and the rest-frame U — V colors as a function of 
stellar mass in Figure 1(b) and (c), respectively. There 
is a tight correlation between K s magnitudes and stellar 
masses. Later, we will limit the sample with M s > 10 10 
Mq considering the i^ s -band morphological analysis limit 
of K S (AB) ~ 22.5 mag (see Section 3 for the detail). The 
galaxies in the stellar mass range are located in redder 
part of the bimodal distribution of the rest-frame U — V 
color in this redshift range. At z = 0.8-1.2, - 80% (60%) 
of the galaxies in the MODS with M s > 3 x 10 10 M 
(M s = 1-3 x 10 10 Mq) and with no X-ray-detection (see 
Section 2.3 for identification of X-ray sources) have spec- 
troscopic redshifts. 

2.2. MIPS 24^rn flux and IR luminosity of the z ~ 1 
Galaxies 

In order to estimate IR luminosities, we looked for 
a Spitzer /MIPS (Multiband Imaging Photometer for 
Spitzer) 24 /im source counterpart for each z ~ 1 galaxy. 
The MIPS 24 /im source catalog provided by the GOODS 
team (M. Dickinson et al. in preparation; R. Chary ct al. 
in preparation) contains sources with a 24 /im flux density 
(JW) larger than 80 /iJy with 80% completeness. The cat- 
alog provides not only coordinates of the 24 /im sources 
themselves but also those of the IRAC counterparts if ex- 
ists. We used the coordinates of the IRAC counterparts 
if available for the search. We identified 48 24 /im coun- 
terparts in total for 808 galaxies with .Kg (Vega) < 24 mag 
and 0.8 < z < 1.2. In Figure 1, the 24 /im-detected sources 
are marked with large filled circles (red). Most of them 
have K S (AB) < 22 mag, M s > 10 10 Mq, and bluer col- 



ors in the rest-frame U — V than the 24 /xm-undetected 
galaxies in the same stellar mass range. 

Among the i^ s -selectcd galaxies with no 24 /jm detec- 
tion, we examine 112 galaxies with M s > 1 x 10 10 Mq (in 
our final sample defined in Section 4.1) on the MIPS 24 
/im image (version 0.3) visually in order to estimate how 
many galaxies would be missed in the 24 /xm catalog by 
confusion effect. We find that 12 out of 112 galaxies are 
located near bright 24 /jm sources, and two of them have 
M s > 3 x 10 10 Mq. If those 12 galaxies have a 24 /im flux 
larger than 80 /iJy, the fraction of 24 /im-detected galax- 
ies in our final sample will increase about 3% (13%) at M s 
> 3 x 10 10 (= 1-3 xlO 10 ) Mq. 

For the 24 /im-detected galaxies, we estimate their bolo- 
metric IR luminosities Lir, (defined as Ltr = L[8-1000 
/mi]) from the observed 24 /im fluxes. Local IR- luminous 
(Xir > 10 10 Lq) galaxies show a tight correlation between 
12 /xm flux and Ltr. estimated from 12, 25, 60 and 100 /xm 
MIR flux densities with a scatter of ~ 0.15 dex (Chary & 
Elbaz 2001). Following their prescription, we convert 
the 24 /im fluxes (in /iJy) of our z ~ 1 galaxies to the 
rest-frame 12 /im luminosities, and then compute their 
bolometric IR luminosities in solar luminosity as follows: 

v e S Va (\2 /mi) = v Q f 2i x AttD l (z) 2 , (1) 
L m = 0.89l°i? x [v e S Ue (U /im)] 1 ' 094 , (2) 

where Dl is the luminosity distance corresponding to the 
redshift z, v Q is the observed frequency corresponding to 
the wavelength A = 24 /im, and v e is the rest-frame fre- 
quency corresponding to the redshift z. The symbol S v is 
the monochromatic flux at a frequency v expressed in erg 
s _1 Hz -1 . We need to note that the IR luminosity esti- 
mated using the rest-frame 12 /im flux density has an un- 
certainty of ~ 0.3 dex due to the uncertainty of the shapes 
of the IR SED (see Lc Floc'h ct al. 2005; Marcillac et al. 
2006). 
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The 24 /jm flux density limit of /24 = 80 fiJy corre- 
sponds to the IR luminosity Lir of 0.6 x 10 11 , 1.1 x 10 11 
and 2.1 x 10 11 L Q , and SFR of - 10, 15, and 35 M yr" 1 
at z = 0.8, 1.0 and 1.2, respectively In this paper, we 
regard all the 24 /zm-detected (/ 2 4 > 80 /^Jy) galaxies as 
LIRGs while 24 ^m-undetected (/ 2 4 < 80 /^Jy) galaxies 
as non-LIRGs. We need to note that the LIRG identifica- 
tion is incomplete at z = 1.0-1.2 in the strict definition of 
LIRGs with L m = 10 n -10 12 L Q . 

2.3. Identifying Chandra X-ray Sources 

A part of the z ~ 1 galaxies is also detected in X-ray; 18 
out of 808 galaxies with K s < 24 mag and 0.8 <z < 1.2 are 
detected in X-ray source catalog from the Chandra X-ray 
Observatory deep survey image (Alexander et al. 2003). 
In Figure 1, those X-ray-detected sources are marked with 
open circles. We checked those sources in the ACS Zsso- 
band image, and all but one source are extended. 

The evaluation of morphologies of those galaxies may 
be affected by bright point source of type I AGN. Also, if 
the SED of a galaxy is strongly contaminated by an AGN, 
no reliable stellar mass could be obtained without consid- 
ering the effect of AGNs in the SED fitting. Moreover, 
MIR fluxes could not be a good indicator of SFR if their 
MIR radiations are mainly powered by UV radiation from 
AGNs. We will indicate those X-ray-detected galaxies 
with different symbols to the others in later figures, and 
exclude from the statistics. 

3. Morphological Analysis 

3.1. Ser sic profile fitting: Method 

We evaluate the rest-frame NIR morphologies of the 
z ~ 1 galaxies using the MODS _RT s -band image which is 
constructed from good seeing images. The MODS deep 
-Ks-band data consists of images with various seeing sizes, 
FWHMs of from 0."4 to If 2. The image used in K09 
was the deepest one created by combining images with 
the FWHM smaller than 0"8, which has a total integra- 
tion time of ~ 28 hours, 3ct limiting magnitude of K S (AB) 
~ 26.4 mag for point sources, and FWHM of the PSF 
of ~ 0."46. Here, for the morphological analysis, we use 
a shallower but sharper image constructed from images 
with the FWHM < 0."5, which gives a total integration 
time of ~ 16 hours and 3er limiting magnitude of i^ s (AB) 
~ 26.0 mag for a point source, and the final FWHM of 
0"40, which corresponds to a physical scale of ~ 3.2 kpc 
at z = 1 in the cosmology we use in this paper. 

In order to describe morphologies of the z ~ 1 galaxies, 
we use a single component Sersic model (Sersic 1968; see 
also Graham & Driver 2005 for detail of the model). It 
describes the radial surface brightness profile of a galaxy 
by a function given by, 



I(r) = J(r )cxp< -b n 



l/n 
- I -1 



(3) 



y-y c 



1/2 



(4) 



where x and y are aligned with the scmimajor and scmimi- 
nor axes, r is the elliptical radial distance at a pixel (x,y) 
from the center of a source (x c ,y c ), q is the axial ratio of 
the semiminor to the semimajor axis radius, r c is the ma- 
jor axis half-light radius or effective radius which contains 
the half of the total flux, and n is the Sersic index which 
determines how the light profile concentrates around the 
center and how the profile extends to the outskirt. The 
b n is a normalization constant which is a function of n 
and is chosen so that r e is equivalent to the half-light ra- 
dius. The Sersic model includes a wide range of profiles 
of local galaxies such as the exponential (n = 1) and de 
Vaucouleurs (n = 4) models. In fact, a tight correlation 
between n and Hubble T-type (T) is observed in the local 
universe (Ravindranath et al. 2004) and even at z ^ 1 
(Pannella et al. 2006); on average, galaxies with n < 
2-2.5 are mostly disk-like objects (T > 2) while galaxies 
with n > 2-2.5 are mostly spheroids (elliptical and bulge- 
dominated galaxies; T < 2). 

We use the two-dimensional surface-brightness profile 
fitting code, Galfit ver.2.0.3c (Peng et al. 2002). 
Galfit convolves a two-dimensional model profile with 
a point-spread function (PSF) defined by user, and min- 
imizes x 2 residuals between the model profile and a 
real galaxy profile with a Lcvcnbcrg-Marquardt algo- 
rithm. In the x 2 minimization, there are 7 free pa- 
rameters: Sersic parameters (n and r ), total magni- 
tude, semiminor-to-semimajor axial ratio q, position angle 
PA, and the center position (x c ,y c ) of a galaxy. Initial 
guesses for those parameters except for n (being set to 
1.5) were estimated based on the output parameters from 
SExtractor (Bertin & Arnouts 1996) object detection 
code. To estimate x 2 correctly, we take into account the 
pixel-to-pixel signal-to-noise ratio. The dominant source 
of noise for NIR imaging data is Poissonian noise from sky 
background radiation. Therefore we use the square root 
of the exposure map as pixel-to-pixel weights in the x 2 
minimization. 

The PSF used in this paper is obtained by stacking 
about 10 spectroscopically-idcntificd isolated stars. The 
stars are distributed over the if s -band image, and their 
FWHMs are quite uniform (3.42 ± 0.02 pixel). We there- 
fore ignored the dependence on the position in the image, 
and combined all the stars to create a composite PSF, 
which is applied to all the z ~ 1 galaxies. 

For a fitting region, one must have a sufficiently wide 
image size of a galaxy to cover the outskirts of the profile 
in the calculation. Following the galaxy two-dimensional 
profile fitting procedure used in the Galaxy Evolution 
from Morphology and SEDs (GEMS) survey (Haussler 
et al. 2007), we use an image stamp for each galaxy 
with a size (height and width) of 2.5 times the semimajor 
axial length of the Kron elliptical aperture derived from 
the SExtractor output which is expected to contain 
more than 90% of the total flux of the galaxy regardless 
of brightness (Kron 1980). Neighboring galaxies around 
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the galaxy of interest were masked and excluded from the 
fitting, or were fitted simultaneously if they are so close to 
the galaxy of interest. We determined whether a neighbor 
should be masked out or fitted simultaneously based on 
their Kron elliptical apertures same as in Haussler et al. 
(2007). First, the apertures are enlarged by increasing 
the semiminor and the semimajor axis radii by a factor 
of 1.5, and then if the neighbor's extended aperture over- 
laps the extended aperture of the galaxy of interest, the 
neighbor is fitted simultaneously; otherwise, the neigh- 
bor is masked out with the extended aperture. Neighbors 
with _R' S (AB) > 22.5 mag are masked out regardless of the 
overlap because the profile fitting for them would not be 
reliable (see the next subsection). Note that a difference 
in the treatment of neighbors (masked or fitted) has little 
effect on the fitting result for a galaxy of interest, except 
for galaxies with very bright neighbors, which suffer flux 
contamination from the neighbors and their fittings fail or 
derive obviously unreliable results with bright residuals. 

3.2. Sersic profile fitting: Simulation 

To quantify how reliably we evaluate morphology of 
galaxies with the MODS i^ s -band image, we perform a 
Monte Carlo simulation with ~ 1000 artificial galaxies. 
Galfit is used to generate artificial galaxy images with 
parameters given randomly in the following ranges: 18 S. 

^ 25 pixel (equivalent to 
: = 1), 0.5 < n(K s ) < 10, 



< 



K S (AB) < 24 mag, 0.3 
0."035-2"9 or 0.3-23 kpc at 

0.1 ^ q £ 0.9, and 5j PA £ 180 deg. Each artificial 
galaxy is convolved with the PSF used for the real galax- 
ies, and is embedded randomly in the if s -band image. 
S Extractor is used to detect them and measure appar- 
ent shapes, which are used as initial guesses for Galfit 
in the same manner as for the real galaxies. 

Using the number of galaxies detected among galax- 
ies embedded, we first show detection completeness as a 
function of input magnitude and effective radius for low- 
ninput(-Ks) (n input [K s }< 2.5) sample and high-?i input (i ; Q 
("input [K s ] > 2.5) sample in Figure 2(a) and (b), respec- 
tively. Also shown in both panels is the distribution of the 
K s magnitude (MAG_AUTO from SExtractor) and ef- 
fective radius evaluated with Galfit of the z ~ 1 galaxies 
with M s > 1 x 10 10 M Q according to n(K s ). Although 
faint and large galaxies (located at top-right in Figure 2) 
can be missed at K S (AB) > 22.5 mag in the low-n(K s ) 
sample, all the detected z ~ 1 galaxies are located well 
below the detection limit. 

Comparing input parameters with measured ones for 
the detected artificial galaxies, we estimate the uncer- 
tainty of the profile fitting parameters. For this analysis, 
we only select the artificial galaxies having measured K s - 
band magnitudes and effective radii within 20% of those 
of the observed z ~ 1 galaxies (shown in Figure 2 with 
filled symbols). Figure 3 shows the differences between 
the input and measured parameters as a function of K s 
input magnitude. There is no systematic offset in any pa- 
rameters up to -Ks. i n p U t (AB) =22.5 mag. At i^ s .i nput (AB) 
~ 22.5 mag, we find that low-ni npu t(if s ) galaxies tend to 
have — 1(±34)% lower Sersic index, 0(±9)% lower size, 
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Fig. 2. Detection completeness map for galaxies in the 
MODS ii's-band image estimated from the artificial galaxy 
simulation, (a) The map for artificial galaxies with n nlpu t(Ks) 
< 2.5. The grey-scale map indicates the simulation result. 
Ovcrplottcd large (small) circles show the observed effec- 
tive radii and K B magnitudes of the z ~ 1 galaxies with M B 
> 3 X 10 10 Mq (M s = 1-3 Xl0 10 M Q ) and n(K B ) < 2.5 in 
the MODS, (b) Similar to the panel (a), but for galaxies with 
n(K B ) > 2.5. 
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Fig. 3. Magnitude dependences of the difference between in- 
put and output structural parameters for artificial galaxies 
(circles for those with rainputl-^s] < 2.5 and triangles for those 
with fiinputl^s] > 2.5). From the top to bottom panel, the 
median difference and Icr scatter are shown for Sersic indices 
n(K B ), effective radii r e (K B ), and axial ratios q(K B ). Among 
the galaxies detected, we used only those with the similar 
K B magnitude (±0.2 mag) and the effective radius (80-120%) 
compared to the observed z ~ 1 galaxies. 
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Fig. 4. Difference between input and output values for each 
structural parameter as a function of input value for artificial 
galaxies. From the top to bottom panel, the median difference 
and la scatter are shown for Sersic indices n(K s ), effective 
radii r e (K B ), and axial ratios q(K B ). The symbols are the 
same as in Figure 3. Among the galaxies used in Figure 3, we 
used only those with i^ s ,input(AB) < 22.5 mag. 



and 1 (±14)% larger axial ratio, and higIi-ni, lpu t(.K s ) galax- 
ies tend to have — 3(±45)% lower Sersic index, — 4(±29)% 
lower size, and — 5(±16)% lower axial ratio. 

Finally, we examine differences between input and out- 
put values for each parameter as a function of input value 
of the galaxies with K S (AB) < 22.5 mag, which are shown 
in Figure 4. We find no significant trend in each pa- 
rameter for the low-ni n p U t(if s ) sample. On the contrary, 
for the high-ni npu t(.Ks) sample we find a systematic off- 
set of An/n ~ —15% and Ar c /r c ~ —10% for galaxies 
with ni npu t(.Ks) ^ 8. Although those offsets less affect 
to our morphological separation criterion of n(K s ) =2.5 
(described in Section 4.1), we have to keep in mind the 
offset in size when we study the sizes of the high-n(K s ) 
galaxies. 

Based on the results of the simulation with artificial 
galaxies, we find that the morphological analysis using 
Galfit with the MODS deep ^C s -band image is reliable 
up to K S (AB) ~ 22.5 mag. Therefore, in later sections we 
will limit our z ~ 1 morphological sample at the stellar 
mass of M s ~ 10 10 Mq corresponding to that magnitude. 

3.3. Non-parametric morphological analysis 

We perform an another morphological analysis based 
on a non-parametric method: Concentration- Asymmetry 
(C-A) analysis (Abraham et al. 1994; Abraham et al. 
1996; Conselice et al. 2000; Conselice 2003). Unlike 
Galfit, the C-A analysis does not require information 
on the shape of the PSF and assumptions on the shape of 
the light profile. Therefore, we use this method in parallel 
with the Sersic profile fitting. 

The concentration index (C) is an indicator how the 
light profile of a galaxy concentrates around the center. 
The definition used in this paper is, 



C = 51og! 



rgo 



(5) 



where r^o and T2o are the radii which contain 80% and 
20% of the total flux, respectively. The total flux is defined 
as MAG_AUTO derived from the SExtractor output. 
This parameter is related to n of Sersic profile and there 
is certainly a correlation between these two parameters. 

The asymmetry index (A) is an indicator how the two- 
dimensional light profile is disturbed, which is defined as, 



S| B — Bisq\ 
Sliol 



,(6) 



where Jo and Jiso are the intensity of each pixel of the 
image and of the image rotated by 180° around the galaxy 
centroid, respectively. The symbols Bq and i?iso have the 
similar definitions but for the background region. The A is 
defined as the minimum value with varying the position of 
the center of rotation around the original galaxy centroid. 
The summation taken in the area is defined using the Kron 
elliptical aperture used in Section 3.1. Wc found that the 
contribution of sky background, the term ~S\Bq — BisoL 
correlates well to the elliptical aperture size, and is less 
dependent on the position on the image. Therefore, we 
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Fig. 5. K s magnitude distribution of z ~ 1 galaxies as a func- 
tion of stellar mass. The filled symbols show the LIRGs while 
open symbols show the non-LIRGs. The sample is divided by 
the Sersic index: circles for lovf-n(K B ) (n[K B ] < 2.5) galax- 
ies and triangles for high-n(A" s ) (n[K s ] > 2.5) galaxies. The 
X-ray-detected sample is marked with black dots. The hori- 
zontal dashed line represents K B (AB) = 22.5 mag, which in- 
dicates the limit of morphological analysis. This magnitude 
corresponds to the stellar mass of ~ 1 X 10 10 Mq which is 
adopted as our stellar mass-limited morphological sample cri- 
terion. The vertical dashed line separates less massive (M s = 
1-3 xlO 10 Mq) and massive (M s > 3 X 10 10 Mq) galaxies. 
The less massive galaxies are indicated with smaller symbols. 

assume a universal background over the image and scale 
it according to the aperture size of galaxies. 

4. Results 

4-.1. Rest-frame NIR morphologies of z ~ 1 LIRGs 

The morphological simulation shows that the MODS 
-ftT s -band image allows us to evaluate a light profile of 
galaxies with X S (AB) ~ 22.5 mag or brighter (Section 
3.2). As seen in Figure 1(b), this K s magnitude limit 
corresponds to M s ~ 10 10 Mq at z ~ 1. We therefore ap- 
ply a stellar mass cut of M s = 1 x 10 10 Mq to the z ~ 1 
galaxy sample to define a morphological sample, which 
contains 139 galaxies in total. Considering the scatter in 
K s magnitude at M s = 1 x 10 10 Mq (representing a vari- 
ety of A/ s /L), we also construct a more conservative mor- 
phological sample containing galaxies with M s > 3 x 10 10 
M Q , to which 67 out of 139 galaxies belongs. We refer 
the galaxies with M s > 3 x 10 10 M Q (M s = 1-3 xlO 10 
Mq) as massive (less massive) galaxies. Our morphologi- 
cal sample achieves high completeness thanks to the very 
deep imaging data; 100% and 96% for the massive and 
less massive sample, respectively. We summarize in Table 
1 the numbers of the z ~ 1 galaxy morphological sub- 
sample in the MODS deep region for stellar mass cuts 



Fig. 6. Fraction of LIRGs (galaxies with /24 > 80 /^Jy) as a 
function of stellar mass. The samples are split into two groups 
(n[.Ks] < 2.5: circles and n[K B ] > 2.5: triangles). All the 
X-ray-detccted galaxies are excluded. The error bars account 
for Poissonian errors only. 

of M s = 1-3 xlO 10 Mq and M s > 3 x 10 10 Mq. In the 
parentheses, the numbers of galaxies with spectroscopic 
redshift are shown. The spectroscopic redshifts are avail- 
able for about 60% and 80% of less massive and massive 
galaxies, respectively. The uncertainty of the photomet- 
ric redshift estimates of 0.072 (Section 2.1) brings about 
1-2% uncertainties in physical size (R e ) estimates, which 
is sufficiently small for our study. In Appendices we show 
the properties and images used in this study for all the 
z ~ 1 sample. 

Figure 5 summarizes the morphological evaluation for 
the z r~j 1 galaxies in the stellar mass-its magnitude di- 
agram. The sample is divided by their IR luminosities 
(LIRGs and non-LIRGs) and Sersic indices (less concen- 
trated system: low-n[-ff s ] and highly concentrated system: 
high-n[-K's]) in the figure. The X-ray-detected sample is 
marked with black dots. For M s < 3 x 10 10 Mq, the sam- 
ple is shown with small symbols. In addition to the fact 
that the LIRGs are popular at M s > 10 10 Mq and bright in 
the if s -band as already seen in Figure 1(b), most of them 
are classified in n(K s ) < 2.5. In Figure 6, these trends are 
shown quantitatively. Although with small sample size of 
the LIRGs at the most massive bin, the fraction of LIRGs 
in the sample increases with stellar mass, and most of 
them have n(K B ) < 2.5 while the LIRGs with n(K s ) > 2.5 
are rare. 

Figure 7(a) and (b) show the distribution of the struc- 
tural parameters, the Sersic indices n{K s ) and the physi- 
cal sizes (effective radii in kpc) along the semimajor axis 
R e (K s ), for galaxies with M B = 1-3 xlO 10 M Q and with 
M s > 3 x 10 10 Mq, respectively. In each figure, the bot- 
tom and right panels show the respective distributions 
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Table 1. The number of sample galaxies in this paper* 





M s = 1-3 x 10 iu M Q 
No X-ray All 


M s > 3 x 10 1U M Q 
No X-ray All 


LIRGs with n[K s }< 2.5 
LIRGs with n[K s ]> 2.5 
Non-LIRGs with n[K s }< 2.5 
Non-LIRGs with n[K s }> 2.5 


7 (6) 7 (6) 
2 (1) 3 (2) 
43 (29) 43 (29) 
20 (6) 22 (8) 


27 (21) 30 (24) 
1 (1) 3 (3) 
15 (8) 18 (11) 
24 (21) 29 (26) 


Total 


72 (42) 75 (45) 


67 (51) 80 (64) 



*In the parentheses, numbers of galaxies with spectroscopic redshift are given. 
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Fig. 7. Structural parameters (Sersic indices and effective radii in physical scale along the semimajor axis) measured in the MODS 
deep i\s-band image for the z ~ 1 galaxies with (a) M B = 1-3 xlO 10 Mq, and with (b) M B > 3 X 10 10 Mq. The symbols are the same 
as in Figure 5. The error bars indicate the uncertainties of the fitting. The vertical dashed line separates less concentrated and highly 
concentrated systems (n[ifg] < 2.5 is considered as a less concentrated system). In each figure, the bottom and right panels show 
numbers of the samples at a given stellar mass range as a function of n(K B ) and R C (K B ), respectively. Thick (thin) line corresponds 
to the numbers of the LIRGs (non-LIRGs), and solid (dotted) line corresponds to the numbers of the low-n(K B ) (high-n[Ks]) sample. 
The error bars of the histograms account for Poissonian errors only. All the X-ray-detected galaxies are excluded from the counts. 



of the parameters, where the X-ray-detected galaxies are 
excluded from the counts. Most of the LIRGs are well 
described with exponential disk- like profiles ({n[K s ]) ~ 
1.2). With our separation criterion of n(K s ) — 2.5 in the 
massive (M s > 3 x 10 10 M@) sample, 96% (27/ [27+1]) of 
the LIRGs are classified in less concentrated (\ow-n[K s ]) 
system while only one galaxy (~3%) is in highly concen- 
trated (high-n[if s ]), or bulge-dominated system. As seen 
in the bottom panel of Figure 7(b), the LIRGs have a 
skewed n(if s )-distribution where most of the LIRGs are 
distributed between n(K s ) = 1-2, and therefore the frac- 
tion is not changed significantly even if we use a different 
separation criterion; about 89% and 100% for the criterion 
n(K s ) =2.0 and 3.0, respectively. Melbourne et al. (2008) 
reported the similar result to ours by AO-supported NIR 
observations that two thirds (~ 67%) of z ~ 1 LIRGs are 
disk-like galaxies according to their visual morphological 
classifications. Despite various differences between our 



study and Melbourne et al. (2008) such as the spatial 
resolution of data (seeing-limitcd or AO-supported) or the 
morphological classification methods (based on the profile 
fitting or visual appearance), our study confirms quanti- 
tatively with larger number of sample that most of the 
z ~ 1 LIRGs have disk-like morphologies. 

If we examine the fraction among the massive galaxies, 
we find that about 64% (27/[27+15]) of the \aw-n{K a ) 
galaxies are LIRGs among the galaxies with M s > 3 x 10 10 
Mq . As seen in Figure 6 as well, the LIRG fraction in the 
low-ri^-Kg) sample increases with stellar mass. At z ~ 1, 
the disk-like massive galaxies have surprisingly high frac- 
tion of LIRGs. In fact, the residual images after subtract- 
ing best fit Sersic component show clear spiral structure 
for some of z ~ 1 galaxies. 

As for the sizes, the right panel of Figure 7(b) shows 
that the LIRGs have larger R e (K s ) than the non-LIRGs 
on average. Comparing the \ow-n(K s ) galaxies in between 
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LIRGs and non-LIRGs, a mean size and a sample scat- 
ter of the LIRGs is 3.8 ± 1.6 kpc while that of the non- 
LIRGs is 2.6 ± 1.0 kpc. The high-n(A" s ) non-LIRGs have 
a wide i? (A s )-distribution and some of them have the 
sizes around 1 kpc. In the MODS deep data, a FWHM 
of the PSF is 0"40, which corresponds to R C (K S ) ~ 2 kpc 
at z = 1. Since Galfit takes into account the PSF in the 
profile fitting, sizes smaller than the PSF size can be eval- 
uated in principle. However, the morphological analysis 
of those unresolved galaxies can be easily affected by what 
PSF (size and shape) is assumed in the fitting, and their 
resultant sizes would have larger uncertainties than the 
other larger galaxies. Therefore, we have to pay attention 
to those small (R e [K s ] < 2 kpc) galaxies when discussing 
their sizes. However, for the LIRGs, most of them are 
much larger than the PSF size and their profiles must be 
derived reliably 

4-2. Rest-frame NIR Stellar Mass-Size Relation 

We show the stellar mass-size distribution in the K s - 
band (rest-frame J-band) in Figure 8(a) and (b) for the 
low-n(A s ) and high-n(A s ) sample, respectively. The sizes 



hi 



are shown as circularized effective radii, i? e ,c = Rc x q 



1/2 

where q is the intrinsic semiminor-to-semimajor axial ra- 
tio of the galaxy. The circularized method yields large dif- 
ference in size between the LIRGs and non-LIRGs when 
compared to those derived from the semimajor axial ra- 
dius of the galaxies which is shown in Figure 7. To com- 
pare those sizes with those of local galaxies, we construct a 
local galaxy catalog including spectroscopic redshifts, stel- 
lar masses, and Sersic parameters using the NYU Value- 
Added Galaxy Catalog (VAGC; Blanton et al. 2005), 
which is based on the seventh release (DR7) of Sloan 
Digital Sky Survey (SDSS). The stellar masses are esti- 
mated from the fits to the broad-band (u, g, r, i, and 
z centered at 3540, 4770, 6230, 7630, and 9130 A) pho- 
tometry (Blanton & Roweis 2007) using a Chabrier IMF 
(Chabrier 2003) that we convert to the Salpeter IMF 
by adding a factor of 0.25 dcx to them. The Sersic pa- 
rameters are measured in those five bands with a single 
component Sersic model fitting (Blanton ct al. 2003). 
Considering completeness of the SDSS data (Shen et al. 
2003) and possible selection biases (Franx et al. 2008), 
we use the VAGC sample between redshifts of 0.05 and 
0.07, the r-band apparent magnitudes of 15.0 and 17.7 
mag, and the r-band surface brightnesses brighter than 
23 mag arcsec -2 , which contains a sample of ~ 4 x 10 4 
galaxies in total. The criteria except for the redshift range 
affect weakly to the selection. In practice, we have ana- 
lyzed a sample selected with only the redshift criterion 
(0.05 < z < 0.07), and found negligible change in the size 
distribution of the sample. The median redshift of the 
SDSS galaxies is 0.062 at which the SDSS z-band corre- 
sponds to the rest-frame ~ 8600 A. We overplot in Figure 
8 the effective radii measured in the SDSS z-band of the 
local galaxies as a function of stellar mass. The late- and 
early-type galaxies are separated at the Sersic index in the 
z-band, n(z), = 2.5. We see for the both morphological 
types the LIRGs are comparable in size to or slightly (~ 



n(K s ) < 2.5 



-O 1> 



0.0 



0.2 



0.4 0.6 
Axial Ratio q(K s ) 



0.8 



1.0 



Fig. 9. The histogram of axial ratio measured in the 
-ft's-band for the z ~ 1 galaxies with n(K B ) < 2.5 in the 
MODS. The symbols are the same as in Figure 5. Only the 
galaxies with M s = 3-15 xlO 10 Mq are included, and the 
X-ray-detected galaxies are excluded. The error bars account 
for Poissonian errors only. 

20% at a maximum) smaller than the local galaxies while 
the non-LIRGs are significantly (30-70%) smaller than 
both the z ~ 1 LIRGs and the local galaxies. In particu- 
lar, about a half of the low-n(A s ) LIRGs have comparable 
sizes to the local disk-like galaxies. In the local universe, 
a size dependence on the measured wavelength, namely a 
color gradient, between 11000 A and 8600 A is known to 
be only a few per cent (smaller in redder band) for both 
disk-like galaxies (de Jong 1996; Barden et al. 2005) 
and early-type galaxies (Mcintosh et al. 2005). Even if 
we take into account this color gradient to estimate sizes 
of the local galaxies at A = 11000 A, our result would be 
unchanged. 

The trend that the non-LIRGs are further apart from 
the local relation compared to the LIRGs is in good agree- 
ment with a study by Perez-Gonzalez et al. (2008b), who 
used deeper MIPS 24 fj,m image than ours and found that 
24^m-undetected (/24 < 15 /^Jy) sources are smaller than 
other 24 ^m sources for galaxies with M s > 10 11 M . We 
find that the larger difference in circularized size between 
the low-n(A s ) LIRGs and low-n(A s ) non-LIRGs is partly 
caused by the fact that those two groups have different dis- 
tributions of axial ratio, which is shown in Figure 9. For a 
comparison of axial ratio, here we limit the sample within 
AI S = 3-15 xlO 10 M Q which is the stellar mass range of 
the low-n(AT s ) non-LIRGs. A median axial ratio and a 
scatter for the LIRGs and non-LIRGs arc 0.57 ±0.21 and 
0.38 ±0.26, respectively. Taking into account the differ- 
ence in axial ratio, we rcplot the stellar mass-size relation 
at z <~ 1 using not the circularized sizes (i? .c) but the 
semimajor axis sizes (R e ) in Figure 10. The LIRGs are 
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Fig. 8. The circularized effective radius of the z ~ 1 galaxies with (a) n(K B ) < 2.5, and with (b) n(K B ) > 2.5 as a function of stellar 
mass. The symbols are the same as Figure 5. The error bars indicate the uncertainties of the fitting. The median and dispersion of 
the distribution of effective radius in z-band of local late-type (n[z] < 2.5) and early-type (n[z] > 2.5) galaxies (redshifts of ~ 0.062) 
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Fig. 10. Similar to Figure 8(a) but the sizes are described as 
the effective radii measured along to the semimajor axis. The 
symbols are the same as in Figure 5. The error bars indicate 
the uncertainties of the fitting. 



still larger than the non-LIRGs although the size differ- 
ence between the LIRGs and non-LIRGs seen in Figure 8 
becomes small. 

Although some of the law-n(K s ) non-LIRGs with low 
axial ratio (q[K s ] < 0.4) show redder U — V colors than 
the \ow-n(K s ) LIRGs which may be caused by their high 
inclination angles, we cannot investigate further the cause 
of the different q(K s ) distributions considering the small 
sample size. 

4-3. Concentration- A symmetry indices 

Figure 11 shows the result of Concentration- Asymmetry 
(C-A) measurement in the MODS deep i4T s -band image of 
the z ~ 1 galaxies. Since the wavelength in which the 
indices measured is longer than that in previous (optical) 
studies, separation criteria of morphological type of galax- 
ies on the C-A plane established by those studies (e.g., 
Conselicc 2003) are not applicable on our (rest-frame J- 
band) C-A plane. Therefore, it is difficult to interpret 
absolute values of C and A measured in the i^s-band. 
Nonetheless, relative comparison of those values within 
the z ~ 1 sample (LIRGs versus non-LIRGs or \ow-n(K s ) 
versus high-n(A" s ) sample) is meaningful to investigate the 
difference between the populations. Figure 11 shows the 
consistent distribution of morphologies with the result of 
the Galfit analysis. Most of the LIRGs have low C(K S ) 
compared to the non-LIRGs which are distributed widely 
in C(K S ). That is the similar trend as seen in the clas- 
sification with Sersic index. We actually find that C(K S ) 
correlates with n(K s ) well in our sample. 

A part of galaxies show larger A(K S ) than the others. 
At the index A(K S ) > 0.1, there are 12 galaxies (~ 15% 
of the total) for the massive sample. High A(K S ) indices 
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Fig. 11. Asymmetry (A) and Concentration (C) indices of 
the z ~ 1 galaxies measured in the MODS deep its-band im- 
age. The symbols are the same as in Figure 5. The dotted 
lines connect probable physically-interacting pairs in the mas- 
sive (M s > 3 X 10 10 Mq) sample where two galaxies have close 
spectroscopic redshifts as well as being located closely (< 1.5 
arcsec) in the A" a -band image. The galaxy ID:0072 has a com- 
panion (ID:0140) which is in the less massive (M a = 1—3 X 10 10 
Mq) sample. The galaxy ID:2769 (z p hot = 0.98) shows a high 
A due to the existence of a companion (ID:2722) which have a 
spectroscopic redshift of 0.581 (not shown in the figure); This 
pair is just due to a random projection. The galaxy ID:3294 
is located close to very bright source. The galaxy ID:3954 is 
located close to the edge of the image. 

of two out of 12 galaxies are partly affected by nearby 
bright source (ID:3294) or by located near the edge of 
the image (ID:3954). The remaining 10 galaxies have a 
companion very close (< 1.5 arcsec or < 12 kpc in pro- 
jected distances) to themselves. Among the 5 pairs, 4 
pairs have almost the same spectroscopic redshift to each 
other, which means that they can be physically interact- 
ing, not a random projection. Those interacting pairs are 
connected with dotted lines in Figure 11. The remaining 
one pair might be just due to the projection; the galaxy 
ID:2769 have a photometric redshift of 0.98 while the com- 
panion have a spectroscopic redshift of 0.581. Therefore, 
highly asymmetric features observed in the z ~ 1 sample 
are thought to be induced mainly by strong interaction 
(i.e., major-merging events). IR luminosities (an indicator 
of dusty star formation activity) are shown as a function 
of the asymmetry indices in Figure 12. The median and 
l-o" scatter of the index A(K S ) for the LIRGs and the non- 
LIRGs are 0.039 ± 0.076 and 0.034 ± 0.105, respectively. 
We find no correlation between star formation and NIR 
asymmetry property. 



Fig. 12. IR luminosities of the z ~ 1 galaxies in the MODS 
as a function of Asymmetry index measured in the if B -band. 
The symbols are the same as in Figure 5. The non-LIRGs 
(i.e., 24 fim- undetected galaxies) are also marked with down- 
pointing arrows as their upper limits of Ltr corresponding to 
/24 = 80 fi3y at their redshifts. 

5. Discussion 

As seen in Section 2.2, our sample contains majority 
of z ~ 1 LIRGs which are thought to be a major star- 
forming population at that epoch. Quantitative morpho- 
logical analysis reveals that those LIRGs show \ow-n(K s ) 
(disk-like) structure in the NIR. This is consistent with 
the result of the AO-supported study in the if -band of 
Melbourne et al. (2008). Also, the \ow-n(K s ) LIRGs con- 
sist of more than a half of the whole low-n(K s ) galaxies at 
z ~ 1 above M s > 3 xlO 10 Mq. Those observational evi- 
dences indicate that the star formation at z ~ 1 is mainly 
occurred in relatively massive disk-like galaxies, and that 
those star-forming massive disk-like galaxies are popular 
at z <~ 1. In addition, the fact that the \ow-n(K s ) galaxies 
in our sample does not show any correlation between the 
asymmetry and the IR luminosity except for a few close 
pairs indicates that the star formation is not triggered 
mainly by early-stage major-merging events which disturb 
the NIR morphology strongly while it is still possible that 
some of the LIRGs are in late stage of interaction. As 
a next step, we investigate rest-frame UV-to-optical mor- 
phologies of the low-n(-ftTs) galaxies to understand how the 
star formation is occurred in the massive disk-like galax- 
ies. 

We use the ACS Vgo6~ and Zgso-bands (rest-frame U- 
and S-bands) images and the source catalogs version 2.0 
provided by the GOODS project to examine the distri- 
bution of star-forming regions in the z ~ 1 galaxies, and 
to compare them with those of local disk-like galaxies. 
First, we find a zsso-band counterpart for all the z ~ 1 
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Fig. 13. The comparison of effective radii measured in the 
K s - and zgso-band. The symbols are the same as in Figure 5. 
The dotted diagonal indicates where the parameters derived 
in the zsso-band are equal to those in the .ffs-band. The error 
bars indicate the uncertainties of the fitting. 

(-K's-selected) galaxies. For galaxies with multiple zg§Q- 
band counterparts, we select the closest one as a counter- 
part. Two of the 139 z ~ 1 galaxies (ID=3294, 3844) have 
no Zg5o-band counterparts. Both of them are close to a 
bright star, which may have made the Zsso-band source 
detection failed in those region. We do not try to tune 
detection parameters to detect them, but remove them 
from the current discussion. Then, we perform the profile 
fitting for the zsso-band counterparts using Galfit in a 
similar manner as for the i<T s -band galaxies except for fit- 
ting the background sky level simultaneously in the zsso- 
band image. Similar to the if s -band morphological anal- 
ysis, we use the SExtractor outputs as initial guesses 
except for n which is set to 1.5. We use a weight map 
produced by the GOODS team for the signal-to-noise per 
pixel. We create a PSF image for each ACS tile by stack- 
ing spcctroscopically-identificd, unsaturated stars in the 
tile. We confirm that even if we smooth and resample the 
Z850-band image to match the image qualities (i.e., pixel 
scale and PSF size) to those of the if s -band image, the 
fitting results are little changed in a statistical sense in 
comparison with results evaluated with the original zsso - 
band image. Then, we evaluate the V(306-band morpholo- 
gies of the galaxies in the similar manner. Among the 
Z850 _ band structural parameters, the scmimajor, semimi- 
nor axial radius, the axial ratio, and the position angle are 
used as initial guesses for Galfit. The position x, y are 
fixed at those in the zsso-band. The initial guess of the 
total magnitude is obtained from the GOODS-N public 
catalog. 

Figure 13 shows the comparison of the effective radii 
along the scmimajor axis measured in the _ftT s -band and 



zgso-band. We mention that there is no significant differ- 
ence in Sersic index and axial ratio between the rest-frame 
UV-to-optical and NIR bands. The R e (K s ) is systemati- 
cally smaller than the R e (zs5o) irrespective of the IR lu- 
minosity; a mean size and an uncertainty of the mean of 
R e (K s )/R e (z 850 ) ~ 0.65±0.01 and 0.44±0.01 for the low- 
and high-n(K s ) massive sample, respectively. The differ- 
ence between rest-frame J-band and B-baiid sizes found 
in our z ~ 1 sample is a factor of 3-5 larger than that 
for local galaxies. Similar trend of larger color gradient 
of z ~ 1 galaxies was reported by Trujillo ct al. (2007). 
They used the ACS F814W and Near-Infrared Camera 
and Multi-Object Spectrometer (NICMOS) F160W for 
the size comparison of 27 galaxies with M s > 10 11 M Q 
at 0.8 < z < 1.8, and found that sizes measured in F160W 
are 19 ±7% smaller than those in F814W. Previous optical 
studies which investigated sizes of disk-like galaxies with 
M s = 10 10 -10 11 Mq in the rest-frame optical band report 
that there is no size evolution from z ~ 1 to the present 
(Bardenetal. 2005; Trujillo et al. 2007). Our data con- 
firms their result; the circularized sizes of the z ~ 1 low- 
n(K s ) sample measured in the ACS zsso-band (rest-frame 
4250 A) are comparable with those of the local galaxies 
measured in the SDSS <?-band (rest-frame 4360 A) in the 
same stellar mass range. Combined with the systematic 
size difference seen in the rest-frame J-band (Figure 8), 
the z ~ 1 disk-like galaxies seem to have a steeper color 
gradient from the local disk-like galaxies. We show the 
comparison of the color gradient of disk-like galaxies at 
z ~ 1 and z = in Figure 14. For the z ~ 1 galaxies with 
n(K s ) < 2.5, the sizes in the Vqqq- and zgso-band normal- 
ized at the -ftT s -band are plotted. The dotted line is an 
observed color gradient of local disk-like galaxies (Bardcn 
et al. 2005) where the slope, a = —0.184, is obtained 
using the VAGC Sersic parameters in a similar manner to 
us. The intercept of the dotted line is determined to make 
the size at J-band equal to unity. Open squares show the 
relative sizes of the local disk-like (n[z) < 2.5) galaxies in 
our SDSS/VAGC catalog (used in Section 4.2) normalized 
at z-band and shifted vertically to make the size at z-band 
lie on the dotted line. We find a steeper slope of the color 
gradient at z ~ 1 than at z = by a linear fit with the inter- 
cept fixed to unity at the J-band. The slope at z ~ 1 de- 
rived using the ACS VI306, z S50 , and the MOIRCS if s -band 
sizes is a = -0.77 ±0.10 (-0.97 ±0.12) and -0.98 ±0.03 
(-0.56 ±0.03) for the LIRGs and the non-LIRGs in the 
massive (less massive) sample, respectively. On the con- 
trary, the slope at z = derived from u. g, and z-band sizes 
of our VAGC catalog is a = -0.19 ±0.00 for the both stel- 
lar mass ranges. We note that the linear fit with all the 
five SDSS bands yields the same slope of -0.19 as with the 
above three bands. Consequently, the z ~ 1 galaxies with 
M s > 10 10 Mq have 3-5 times steeper color gradient than 
the local galaxies at the similar stellar mass range. The 
change in the color gradient between z ~ 1 and z = can be 
interpreted as that while the distribution of young stellar 
population indicated by the rest-frame UV-optical sizes is 
unchanged from z^ 1 to the present, the structure consist- 
ing of old stellar population indicated by the rest-frame 
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Fig. 14. The normalized effective radii as a function of 
rest-frame wavelength of the z ~ 1 disk-like galaxies. Sizes are 
normalized at J-band (dashed line). The symbols except for 
those with error bars and black dots are the same as in Figure 
5. The large circles with error bar represent the medians and 
scatters. Two stellar mass ranges (Af s > 3 X 10 10 Mq, M s = 
1-3 xlO 10 Mq) are shown separately. All X-ray-detected 
galaxies arc excluded. Galaxies whose normalized radius is 
beyond the vertical axis scale are plotted with an up-pointing 
arrow. The squares with error bar show the median and l-cr 
scatter of the radii for local disk-like galaxies used in Section 
4.2. The dotted line represents the color gradient (a slope 
a = —0.184) for local disk-like galaxies (Barden et al. 2005; 
de Jong 1996) where we set the intercept to make the size at 
J-band equal to unity. 

NIR sizes is still being constructed at z ~ 1. Additionally, 
changes in dust distribution in a galaxy should also con- 
tribute to the decrease of the color gradient. 

Recent NIR Integral Field Unit (IFU) observations in- 
vestigating the kinematics of distant (z ^ 0.6) disk galaxies 
revealed that they show a similar Tully-Fisher relation to 
the local disk galaxies (Flores et al. 2006; Puech et al. 
2008; Lemoine-Busserolle & Lamareil 2010). In addition, 
Ncichel et al. (2008) investigated a radial distribution 
of color and SFR along disk of rotating spiral disk galax- 
ies which show similar morphology and kinematics to the 
local disk galaxies, and found that the very active and re- 
cent star formation is occurred in the outer parts of the 
disks. They concluded the observational evidence as a 
rapid, inside-out disk formation of disk galaxies. Our re- 
sults are consistent with such a disk formation scenario, 
in which local disk galaxies form their disk structure from 
the center to the outskirt. Intense star formation activ- 
ity observed as LIRGs at z ~ 1 may be mostly related to 
the formation and growth of disk structure in relatively 
massive galaxies, and contribute significantly to the high 
cosmic SFR density at z ~ 1. 



6. Summary 

In this paper, we studied rest-frame NIR morphologies 
of galaxies at z = 0.8-1.2 from the deep NIR imaging sur- 
vey (MODS; K09) catalog in a part of the GOODS-N field. 
The AVband image, which is the key data for this work, 
covers ~ 28 arcmin 2 and has a PSF ~ 0"4 (FWHM) which 
corresponds to 3.2 kpc at z ~ 1. Cross-correlating our its- 
band galaxy catalog with the Spitzer /MIPS 24 and 
the Chandra X-ray source catalogs, we identified LIRGs 
and AGN candidates. Using Galfit, we performed a two- 
dimensional light profile fitting of the z ~ 1 galaxies in the 
A s -band with a single component Scrsic model. Our sim- 
ulation with artificial galaxy images showed that our mor- 
phological analysis is reliable up to A S (AB) ~ 22.5 mag 
irrespective of morphology. That magnitude limit corre- 
sponds to M s ~lx 10 10 Mq at z ~ 1 and we adopted that 
stellar mass as a selection criterion. As a non-parametric 
morphological study, we measured the Concentration (C) 
and Asymmetry (A) indices of the galaxies in the K s - 
band. 

We investigated NIR morphological properties of the 
139 galaxies having z = 0.8-1.2 and M s > 1 x 10 10 M . 
The main results about those galaxies are summarized as 
follows. 

• - 90% of LIRGs show disk-like (n[K s ] ~ 1-2) light 
profiles in the A s -band. 

• Those low-n(As) LIRGs consist of 60% of the whole 
low-n(As) sample above M s > 3 xlO 10 M©. 

• About a half of the low-n(A' s ) LIRGs are compara- 
ble in size and the others are slightly (~ 20% at a 
maximum) small compared to local disk-like galax- 
ies. 

• No correlation between the NIR asymmetry proper- 
ties and Lir, (star formation activity) is found. 

Those results indicate that the star formation at z ~ 
1 is (i) mainly occurred in most relatively massive disk 
galaxies, and (ii) is not triggered by early-phase strong 
galaxy-galaxy interactions which disturb strongly the NIR 
morphology. 

In order to investigate how the star formation is oc- 
curred in the massive disk-like galaxies, we compared the 
rest-frame J-band morphologies of the z ~ 1 galaxies with 
the rest-frame U- and B-band ones using the HST/ ACS 
Veo6- and Zsso-band image, respectively. There is no sig- 
nificant difference in Sersic index and axial ratio between 
the rest-frame UV-optical and NIR bands. The compari- 
son of the effective radii shows a ~ 30% systematic offset 
where sizes in redder band are smaller. Although this is 
known as "color gradients" of galaxies in the local uni- 
verse, the gradient we found is 3-5 times steeper than in 
the local universe. Since we found the fact that the rest- 
frame optical sizes of the z ~ 1 galaxies are comparable 
to the local galaxies, this steeper color gradient indicates 
that fundamental disk structure in those massive galaxies 
is still being constructed at z ~ 1. Our results indicate 
not only that more than a half of relatively massive disk- 
like galaxies at z ~ 1 are in violent star formation epochs 
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observed as LIRGs but also that most of those LIRGs 
are constructing their disk structure vigorously. The high 
star formation rate density in the universe at z <~ 1 may 
be dominated by star formation in disk region in massive 
galaxies. 
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Appendix 1. Properties of the MODS z ~ 1 
galaxies 

Table 2 is a sample of the information about the proper- 
ties of the z ~ 1 galaxies analyzed in this paper. Columns 
1-6 list the basic properties of the galaxy derived from 
Kajisawa ct al. (2009, K09). Column 1 indicates the 
galaxy identification number. Columns 2 and 3 are the 
Right Ascension (R.A.) and Declination (Dec.) of the 
galaxy for epoch J2000 in degrees. Column 4 lists the ap- 
parent A s -band magnitude in the AB system. Column 5 
is the measured redshift of the galaxy. Column 6 specifies 
whether the redshift was determined spectroscopically (s) 
or photometrically (p). Column 7 is the stellar mass of 
the galaxy in units of 10 10 A/ Q . Column 8 is the MIPS 24 
/im flux in units of fiJy derived from the 24 /an source 
catalog (M. Dickinson et al. in preparation; R. Chary et 
al. in preparation). Possible 24 /im sources which would 



be caused by confusion effect identified by visual inspec- 
tion on the MIPS 24 ^m image are denoted by an aster- 
isk. Column 9 is the IR luminosity in units of lO 11 ^©. 
Columns 10-14 list the morphological properties of the 
galaxy measured in the its-band. Column 10 indicates the 
value of the Sersic index of the fit to the galaxy. Column 
11 is the effective radius along the semimajor axis of the 
galaxy. Column 12 is the axial ratio of the galaxy. And 
finally, Column 13 and 14 are the concentration and asym- 
metry indices of the galaxy. Galaxies listed in the X-ray 
catalog (Alexander et al. 2003) arc denoted as "X-ray" 
in the Comment. The entire catalog for all the 155 galax- 
ies including the X-ray-detected galaxies analyzed in this 
paper is available as an online material. 

Appendix 2. Images of the MODS z ~ 1 galaxies 

Figure 15 is a sample of postage stamp images and K s - 
band surface brightness profiles of the z ~ 1 sample an- 
alyzed in this paper. In the left panels, ACS zsso-band 
(rest- frame B), A s -band (rest-frame J) images, and K s - 
band GALFIT residual image after model subtraction from 
the -fTg-band image are shown from left to right. North is 
to the top and the East is to the left. The size of each 
image is approximately 5.0 x 5.0 arcsec (~ 40 x 40 kpc). 
The residual image is shown with narrower dynamic range 
than that for the A' s -band image to display the residual 
pattern clearly. Simultaneously-fitted neighbors are sub- 
tracted from the residual image as well with their best-fit 
Sersic models. In the right panels, observed (filled cir- 
cles) and best-fit (solid line) surface brightness profiles 
in the A s -band along the semimajor axis are shown in 
the top. The profiles arc derived with elliptical isophotc 
fitting package ellipse in IRAF. The shape (axial ra- 
tio and position angle) and the center of the ellipse are 
fixed with values derived using GALFIT, and the isophote 
radius is changed from 0.0 to 2.5 arcsec. The observed 
profile is derived from the Af s -band image subtracted the 
simultaneously-fitted neighbors with their best-fit models, 
and is plotted until the radius reaches 1.5 times the semi- 
major axial length of the Kron elliptical aperture, which 
is a criterion for choosing neighbors, or 2.5 arcsec at a 
maximum. Error bars represent 1-a scatter of intensity 
data along the ellipse at a given radius. The dot-dashed 
line shows the profile of the PSF used in GALFIT normal- 
ized to match the brightness with the observed point at 
the center. A residual of the observed profile is shown 
in the bottom. Image lists for all the 155 galaxies includ- 
ing the X-ray-detectcd galaxies analyzed in this paper are 
available as online materials. 
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Table 2. Properties of the MODS z ~ 1 galaxies. 



ID 


R.A. 


Dec. 


A' S (AB) 


Redshift 
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Fig. 15. LIRGs with n(K B ) < 2.5 at z ~ 1 in the MODS deep region. The first two galaxies are shown as a sample. Iimage lists for 
all the 155 galaxies including the X-ray-detected galaxies analyzed in this paper are available as online materials. Left panels: ACS 
zg5o-band (rest-frame B), A s -band (rest- frame J) images, and A a -band GALFIT residual image after model subtraction from the 
A s -band image are shown from left to right. North is to the top and the East is to the loft. The size of each image is approximately 
5.0 X 5.0 arcscc (~ 40 X 40 kpc). The residual image is shown with narrower dynamic range than that for the A" a -band image 
to display the residual pattern clearly. Simultaneously-fitted neighbors are subtracted from the residual image as well with their 
best-fit Sersic models. Right panels: Observed (filled circles) and best-fit {solid line) surface brightness profiles in the A s -band 
along the scmimajor axis are shown in the top. The profiles are derived with elliptical isophote fitting package ellipse in IRAF. 
The shape (axial ratio and position angle) and the center of the ellipse are fixed with values derived using GALFIT, and the isophote 
radius is changed from 0.0 to 2.5 arcsec. The observed profile is derived from the A s -band image subtracted the simultaneously-fitted 
neighbors with their best-fit models, and is plotted until the radius reaches 1.5 times the scmimajor axial length of the Kron elliptical 
aperture, which is a criterion for choosing neighbors, or 2.5 arcsec at a maximum. Error bars represent 1-cr scatter of intensity data 
along the ellipse at a given radius. The dot-dashed line shows the profile of the PSF used in GALFIT normalized to match the 
brightness with the observed point at the center. A residual of the observed profile is shown in the bottom. 
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